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Abstract 

Recently hints of lepton flavor non-universality emerged in the BaBar and LHCb 
experiments. In this paper we propose tests of lepton universality in Ut scattering. 
To parametrize the new physics we adopt an effective Lagrangian approach and 
consider the neutrino deep inelastic scattering processes Ur + N —)■ r + X and 
+ N —)■ p, + X where we assume the largest new physics effects are in the r 
sector. We also consider an explicit leptoquark model in our calculations. In order 
to make comparison with the standard model and also in order to cancel out the 
uncertainties of the parton distribution functions, we consider the ratio of total and 
differential cross sections of tan-neutrino to muon-neutrino scattering. We hnd new 
physics effects that can possibly be observed at the proposed Search for Hidden 
Particles (SHIP) experiment at CERN. 
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1 Introduction 


The flavor sector of standard model (SM) has many pnzzles. A key property of 
the SM gange interactions is that they are lepton flavor nniversah Evidence for 
violation of this property would be a clear sign of new physics (NP) beyond the 
SM. In the search for NP, the second and third generation quarks and leptons could 
be special because they are comparatively heavier and are expected to be relatively 
more sensitive to NP. As an example, in certain versions of the two Higgs doublet 
models (2HDM) the couplings of the new Higgs bosons are proportional to the 
masses and so NP effects are more pronounced for the heavier generations. Moreover, 
the constraints on new physics, especially involving the third generation leptons and 
quarks, are somewhat weaker allowing for larger new physics effects. 

Interestingly, there have been some reports of non-universality in the lepton 
sector from experiments. Recently, the BaBar Collaboration with their full data 
sample has reported the following measurements m- 

BiB D+r-pA 

--V=- ^ , ! = 0.440 ± 0.058 ± 0.042 , 

B{B ^ D+i-ue) 

BiB D*+T-Vr) 

—^ = 0.332 ±0.024 ±0.018 , 1 

where £ = e,fi. The SM predictions are R{D) = 0.297 ± 0.017 and R{D*) = 
0.252 ± 0.003 [H [3], which deviate from the BaBar measurements by 2a and 2.7a, 
respectively. (The BaBar Collaboration itself reported a 3.4 ct deviation from SM 
when the two measurements of Eq. ([T]) are taken together.) This measurement of 
lepton flavor non-universality, referred to as the R{Di*i) puzzles, may be providing 
a hint of the new physics (NP) believed to exist beyond the SM. There have been 
numerous analyses examining NP explanations of the R{Di*i) measurements [11[5]. 

In another measurement the LHCb Collaboration recently measured the ratio 
of decay rates for B~^ —)■ K^£^£~ {£ = e, /i) in the dilepton invariant mass-squared 
range 1 GeV^ < < 6 GeV^ |6]. They found 

Rk = ^ = 0-745to:o?4 (stat) ± 0.036 (syst) , (2) 

is a 2.6cr difference from the SM prediction of Rk = 1 ± 0(10“^) [7]. In addition, 
we note that the three-body decay 5° —)■ K* by itself offers a large number of 
observables in the kinematic and angular distributions of the final-state particles, 
and it has been argued that some of these distributions are less affected by hadronic 
uncertainties [8]. Interestingly, the measurement of one of these observables shows 
a deviation from the SM prediction [9]. However, the situation is not clear whether 
this anomaly is truly a first sign of new physics [10] . 

The tan neutrino, was discovered by the DONuT experiment m which 
measured the charged-current (GG) interaction cross section of the tan neutrino. 
The DONuT central-value results for the scattering cross section show deviation 


R{D) ^ 
R{D*) = 


1 





from the standard model predictions by abont 40% bnt with large experimental 
errors; thus, the measurements are consistent with the standard model. The third 
generation lepton has been explored relatively less than the other two generations 
and in particular there has not been much investigation of Vr properties. One of the 
predictions of the Standard Model (SM) is that gauge bosons couple to the three 
generations of leptons universally. A careful test of this prediction is very important 
and observation of non-universality in the interactions of the lepton families would 
be an important discovery. 

In previous publications we considered new physics in scattering for quasi¬ 
exclusive, resonant and DIS scattering [12]. In those papers we were more focused 
on the error in the extraction of neutrino mixing angles in presence of new physics. 
In this paper we focus on observables that may be measured at a Vr scattering 
experiment. There is a proposed Search for Hidden Particles (SHIP) experiment at 
CERN [T3| which is expected to have a large sample of tan neutrinos which could be 
used to probe new physics in scattering. In our previous work we did not include 
new physics tensor interactions which we consider in this work. In this work we 
will be interested at neutrino energies where the DIS component of the scattering 
process is dominant. 

We start with an effective Hamiltonian description of new physics operators. We 
hx the constraints on the couplings from charged current r decays. We consider the 
decays r —)■ ttVt- and r tttiVt- which are well measured. We hnally consider an 
explicit leptoquark model where both scalar and tensor interactions, with relations 
between the couplings of the two interactions, are present. 

The paper is organized in the following manner. In sec. 2 we introduce the 
effective Lagrangian to parametrize the NP operators, describe the formalism of 
the decay process and introduce the relevant observables. In sec.3 we present our 
results and in sec.4 we present our conclusions. We collect some of our equations in 
Appendices (A, B). 


2 Formalism 


In the presence of NP, the effective Hamiltonian for the scattering process Vr + N ^ 
T + X can be written in the form [T4] . 


n 


eff 


^/2 


(1 + 14) [u^^PLd] [h^PLl^l] + 44 [u-i^PRd] 


+SL[uPLd][lPLVi] +Sr [uPRd] [IPlVi] +TL[u(j^^PLd] [la^^PLVi] (3) 


where Gr = 1.1663787(6) x lQ~^GeV~^ is the Fermi coupling constant, Vud is the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element, Pl,r = (1 T are the 

projectors of negative/positive chiralities. We use = z[7^,7i,]/2 and assume the 
neutrino to be always left chiral. To introduce non-universality the NP couplings are 
in general different for different lepton flavors. We assume the NP effect is mainly 
through the r lepton. The effective Hamiltonian involves the quarks of the hrst 
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generations only. It is possible that the qnarks of the other generations will also be 
affected by new physics. We will not assnme any connection between new physics 
for the different generations of qnarks. The SM effective Hamiltonian corresponds 
to 9l = qr = gs = gp = 0. 

The Hamiltonian in the presence of only scalar and tensor operators can be 
written as, 

Ties = + Bs'y5)d /(I - ^ 5)^1 + Tl - 75 )^ - 75 ) 1 /;] , (4) 

where As = Sr + Sl and Br = Sr — Sl with Sl and Sr are the left and right 
handed scalar conplings and Tr is the tensor conpling. 

We will hrst employ a model independent approach and treat the scalar and 
tensor conplings one at a time. Since, in many realistic models both the scalar 
and tensor conplings may be present, we will consider an explicit leptoqnark model 
where both the scalar and tensor conplings are present. 

The Hamiltonian in the presence of only H ± H operators was considered in onr 
previous work [12]. There the effective Hamiltonian was written in terms of a W 
model, which could arise in extensions of the SM [15], as 

£ = 4r/7/V(97-Pi + 9 "-p«)/‘7 + h.c. (5) 


Integrating out the W leads to 

,2 






f'f p I f'fp \ f 


C = 


AGpVfi f 


Bp I 


fl 




^ Ml 


‘■w 


+ h.c., 

[guahtiPivi] + h.c. 


( 6 ) 


Comparing Eq. |6] with Eq. [3] we have the following relations 

f'f i 

-gi g 

f'f I 

-gk g 


Vr = Mknf'fJPt 




(7) 


2.1 Deep Inelastic Neutrino Nucleon Scattering 

In this section we discuss Deep Inelastic Neutrino Nucleon Scattering with the var¬ 
ious types of interactions. 
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2.1.1 Scalar and Tensor Interactions 


In this section, we first present the total and differential cross sections for the deep 
inelastic scattering (DIS) process 


Vj- -\- N —^ T -|- X —^ ^ X^ 


( 8 ) 


with scalar and tensor interactions. The total differential cross section is written 
in terms of contribntions from the standard model, scalar and tensor operators and 
cross terms as follows 

dcTtot _ da^M da^ ^ da^ dasM.ST <^tTs,T 

dxdy dxdy dxdy dxdy dxdy dxdy 

The differential cross section is given in terms of the cross section amplitude as 
follows 

Here, is the four-momentum of the scattered quark, is the target nucleon 

momentum, and ^ is its momentum fraction. /(^) is the the parton distribution 
function (PDF) inside a nucleon and is the incoming neutrino energy. In the 
deep inelastic scattering we calculate the differential cross section with respect to 
the scaling variables which are dehned as follows 


X 

y 


2 z/’ 

V 


( 11 ) 


where x is the Bjorken variable and y is the inelasticity with q being the four- 
momentum transfer of the leptonic probe and 


u = -p. q = M{Ey - Ei). 

The physical regions for x and y are obtained by Albright and Jarlskog dsiin] 

__ < a: < 1 

2 M{E^ — mt) ~ ~ ’ 


( 12 ) 


(13) 


and 


A - B <y < A + B, 


(14) 


where 


A 

B 


1 

2 

1 

2 


m. 


2MEyX 


m. 


2MEyX 






xM\ 

■ 


(15) 

(16) 


4 



















The terms in Eq. [9] are given as 


dasM 

dxdy 


das 

dxdy 
da-r 
dxdy 
dasM,ST 
dxdy 
das,T 
dxdy 

The functions Fi are given as 

F 2 = 

<? 9 

Fs = 2j2U,{tQ^K,„ ( 18 ) 

9.9 


GlME, 

71 




rrif 

4^ 


)F2 


y. 


+{xy{l - -)-y 


mt 


4ME, 


■)Fs 


mt 


2 ME„ 


5 , 


GpME^. 2 t:> 2 \ / \ 7-1 

-{Al + Bl)y{xy + 


47r 
SGlME, 


2ME„ 


''Tl[y{xy + 


TT 


nil 


2ME, 


■)Fi+ 2(1-9/- 


Mxy 


8 +)-^’ - 


ME,, 


0 , 


2GlME^ ^ ^ n y^ ^“e 

-^;—Tl(Bs - As) hvil - j) - yjjpp 


( 17 ) 


where fq and fq are the parton distribution functions inside a nucleon, Vq^q/ is the 
CKM matrix element, and = —q^. 

One can write the differential cross sections above in terms of different variables 
(f, z/) using Eq. [TT]and the transformation [T7] . 


da 

dxdy 


2MEyl' 


da 

dq^dv 


( 19 ) 


5 
























In the new variables, the differential cross sections can be written in the form 
dasM 


dq^du 


Gl 

SttME?, 


+ mj)W, + + ""')) ^2 


1 


M 2 


Gl 


IGtiME?, 


{2ME,q^-u{q^ + mj))W3- 

{Al + + q^)W„ 


2mjEiy 


M 


W, 


r ^2 

'^F rr 2 


7rMR2 


SE ,,!/, 


Tl 2{mi + q^)W, + (SE^^ - {mj + q^) - 


dcTs 
dq^du 
dcJT 
dq^dp 
das^T 

dq^du AnM^El' 
where the (time-reversal invariant) structnre fnnctions are [T7] 

Fi{x) 


ArriiE, 


M 


M 


G% 


,Tl{Bs - As){2E,Mq^ - {mj + q^)u)Ws, 


( 20 ) 


Wiiq^u) = 
W2{q^,v) = 
W,{q\u) = 
W,{q\u) = 


M ’ 

ME2{x) 

V 

ME^{x) 

V 

MF,{x) 


( 21 ) 


We also define some Lorentz invariant variables in terms of the fonr-momenta of 
incoming nentrino {k), target nncleon (p) and the produced charged lepton {k') in 
the laboratory frame 


Q 2 = _g2 = 

W^ = {p + qf. 


( 22 ) 

(23) 


is the magnitude of the momentum transfer and W is the hadronic invariant 
mass. The physical regions of these variables are given by JT2] 

Went <W < Vs -mt, (24) 

in the DIS region with ITcut = 1-4 — 1.6 GeV, and 




where s = {k + pY and 


Ql(W) = ^ 2^(1 ± i3) 


W"- + — mj) 


(25) 


( 26 ) 


with /5 = Aa (1, mj/s, W'^/s) and A(a, 6, c) = — 2{ab -1- 6c -f ca). In the 

lab frame, s = M^ -|- 2 MEi^. 
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2.1.2 Explicit Leptoquark Model 

Here we will discuss an explicit leptoquark model. Many extensions of the SM, 
motivated by a unified description of quarks and leptons, predict the existence of 
new scalar and vector bosons, called leptoquarks, which decay into a quark and 
a lepton. These particles carry non-zero baryon and lepton numbers, color and 
fractional electric charges. The most general dimension four SU{3)c x SU{2)i x 
U{1)y invariant Lagrangian of leptoquarks satisfying baryon and lepton number 
conservation was considered in Ref [IB]. As the tensor operators in the effective 
Lagrangian get contributions only from scalar leptoquarks, we will focus only on 
scalar leptoquarks and consider the case where the leptoquark is a weak doublet 
or a weak singlet. The weak doublet leptoquark, R 2 has the quantum numbers 
(3, 2, —7/6) under SU{3)c x SU{2)l x U{1)y while the singlet leptoquark Si has the 
quantum numbers (3,1, —1/3). 

The interaction Lagrangian that induces contributions to z/,r+A^ t+X process 

is [111120] 

^ 2 *^ = {g^iUiRLjL + g2RQjLi(^2liR) ^ 2 , 

= {9lLMLi(^2LjL + giR,^iRijR) Si, (27) 


where Qi and Lj are the left-handed quark and lepton SU (2 )r doublets respectively, 
while UiR, diR and ijR are the right-handed up, down quark and charged lepton 
SU{2)l singlets. Indices i and j denote the generations of quarks and leptons and 

_ 'j' 

= C'ljj = is a charge-conjugated fermion field. The fermion fields are 

given in the gauge eigenstate basis and one should make the transformation to the 
mass basis. Assuming the quark mixing matrices to be hierarchical, and considering 
only the leading contribution we can ignore the effect of mixing. 

The couplings in Eq. [27|can be constrained from r decays. Because of the doublet 
nature of R 2 there will be additional term like frqq which do not contribute to r 
decays. They will contribute to tau pair production but are much smaller than the 
SM production and hence do not add any new constraints. 

After performing the Fierz transformations, one finds the general Wilson coeffi¬ 
cients at the leptoquark mass scale contributing to the Ur + N ^ r + X process: 


1 

9\l9\r 

91191r' 

2\/2G pVud 

1 

2Ml^ 

'9\l9\r 

2Ml 

9\l9l¥ 


2\/2G pVud 

00 

_ 1 

J 



(28) 


It is clear from Eq. |2H] that the combination from the weak singlet leptoquark and 
the weak doublet can add constructively or destructively to the Wilson’s coefficients 
of the scalar and tensor operators in the effective Hamiltonian. In this section we 
will also consider the possibilities where both the scalar and the tensor operators are 
present and are of similar sizes. In the most general case both the singlet and doublet 
leptoquarks are present and so both the scalar and tensor operators appear in the 
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effective Hamiltonian. As there is limited experimental information, including both 
the singlet and the doublet leptoquarks will allow us more flexibility in htting for the 
Wilson’s coefficients but this will come with the price of less precise predictions for 
the various observables. We can, therefore, consider the simpler cases when only a 
singlet or a doublet leptoquark is present. In these cases, from Eq. [28]the coefficients 
of scalar operators and the tensor operators are related. One can then consider the 
two cases: 

Case, (a): In this case only the weak doublet scalar leptoquark R 2 is present. In 
this case the Wilson’s coefficients are 


1 

allgm' 

2y/2GFVud 

1 

to 

to 

1 _ 

1 

9\191r' 

2\/2G FVud 

- 1 

CA 

00 

_ 1 


( 29 ) 


Case, (b): In this case only the singlet leptoquark is present and the relevant 
Wilson’s coefficients are 


Sl 

Tl 


2y/2GFVud . 

1 

2V2GFVud V 8Ml 


dlWlR 

2 M|^ 

alhln 


(30) 


In the (3,1, —1/3) case, proton decay can occur in presence of an additional lep¬ 
toquark couplings to two quarks. The proton decay only constraints the leptoquark 
mass and the product of this additional coupling with the coupling considered here 
and so we can choose to turn the proton decay constraint as a constraint on the ad¬ 
ditional coupling involving the two quarks. Moreover, the relevant couplings for the 
considered processes involve the third generation and so proton decay constraints do 
not apply in general to these couplings. Also note the proton cannot kinematically 
decay to a r. 

The relations in Eqs. fl29l [30|1 are valid at the leptoquark mass scale. We have 
to run them down to the r mass scale using the scale dependence of the scalar and 
tensor currents at leading logarithm approximation 


Siirnr) 

TL{mr) 


as{mt) 



_as{mr)_ 


_ as{mt) 

' a simt)' 

It 


_as{mr)_ 


_ as{mt) 




“ TL(mLQ) , 


(31) 


where the anomalous dimensions of the scalar and tensor operators are 75 = —6Gf = 
— 8 , 7 t = 2Gf = 8/3 respectively. Further, the beta function, (3q = 11 — 2n//3 

[ 2 TI |23l [23] and nf is the number of active quark flavors. One can use the equations 
above to run down the couplings from a chosen value of jtilq to the tau mass, m,-. 
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In the presence of only one type of leptoquark, singlet or doublet state, one finds 
that the scalar Sl and tensor Tl Wilson coefficients are related to each other at the 
scale of leptoquark mass, S'L(mLQ) = ±4TL(mLQ). 


2.1.3 1/± y4 interactions 

The DIS differential cross section in the presence of V±A operators with respect 
to the variables {x,y) is given in [12]. Here we write it in terms of the momentum 
transfer, using Eq. [19] as follows 


da: 


SM+(V=I=A) 

dq^dv 


Gl 


MMEl 
1 


((|a'p + |6'p) [m] + q^)Wi 


+ 


2 M 

1 


M 2 


(|af + |6f) {AElM - 4E,u - M{mj + q^))W 2 
Re[a'h'*]{2E,Mq^ - + q^))Ws - ^ (|af + \bf) mjE^W^ ) , 


where the dehnitions are 


(32) 


a' = 1 + 7 ^, 

b' = 1 + 7", 

Y = Vl + Vr, 

Y = Vl-Vr. 


(33) 


3 Constraints on NP couplings 

The scalar couplings Sl and Sr can be constrained by the tan decay channel 
T~{ki) —)■ i'r{k 2 ) + 7r~(q'), while the tensor coupling Tl can be constrained by the 
three-body decay channel r(p) —)■ 7r“(pi) + Y{p 2 ) + ^riPs)- In this section we will 
discuss the constraints. 


3.1 r (/ci)z/r(/c 2 ) + TT {q) 


The hadronic current of the bound state tt can be parametrized as 

{0\dY{^ - Y)u\Triq)) = -iY2Yq^, (34) 

where /tt = (92.4 ± 0.1 ± 0.3) MeV [25] is the pion decay constant. The SM decay 
rate is 


SM 









^rfiv • 


(36) 
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Here ^t/tt = 1.0016±0.0014 [26] is the radiative correction. Further, the SM branch¬ 
ing ratio can also be expressed as [27] 

= 0.607Hr(r" ^ z/^e"he) = 10.82 ± 0.02% , (36) 

while the measured Hr(r“ '^~T^T)exp = (10.91 ± 0.07)% [23]. In the presence of a 
scalar state, the decay rate is 


r 


TT 

S 




(37) 


where 

(0| J(H5 - Bs'y^)u\7i{q)) = iV2Um^Bs. (38) 

In order to obtain the scalar hadronic current above, we have multiplied the SM 
hadronic current flM|) by the sum and difference of the quark momenta and used the 
equation of motion - see Appendix (B). The total branching ratio can be written as 
follows 

= BAJm (1 + (rsf ). (39) 

where 



BRl 

R Dtt ’ 

-D-'LSM 


(40) 


with 


r 


TT 

S 


TUr 


(41) 


Note, the interference term of the SM and the scalar NP term vanishes. 

The allowed region of the couplings are given in the contour plot Fig. [T] for 
the measured t~ —)■ within the 2cr level. We consider theoretical uncertainty 

within the Icr level. First, we assume both couplings, are present and take 

the couplings to be real. Next we assume the couplings are complex and take one 
coupling at a time, as shown in Fig. [21 


3.2 t(p) TT (pi) + 7r^(p2) + i^ripa) 

Here we consider two-pion decays of r. The process is 

t(p) -)■ i^riPs) + 7 r~(pi) -h 7 r°(p 2 ). 

The SM and NP amplitudes are 

Msm = ^^^^(7r"7r°|(i7''(l - - 'y^)ur, 


M'j' = 


-iGpVud 


Tl{t^ 7r°|dcr''''(l - - ^^)ur. 


(42) 


(43) 

(44) 
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Figure 1: The constraints on the scalar couplings Sl,r- The colored region is allowed. 
The constraint is from t~ —)■ We treat Sl and Sr as real couplings. 



Re[SJ Re[SR] 


Figure 2: The constraints on the scalar couplings Sl,r- The colored region is allowed. 
The constraint is from t~ —?• Left panel: we take Sr = 0 and treat Sr as a 

complex coupling. Right panel; we take Sr = 0 and treat Sr as a complex coupling. 

We can parametrize the relevant form factors as, 

(7r-7rV7^(l - 7 ^)m|0) = V2F{Q^)kf^, (45) 


{7r-7r°\daf^’'{l - 7'^)n|0) = (46) 
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where k = pi — p 2 and q = pi + P 2 - The form factor F{Q^), along with its error, 
is given in [281 129]. In onr analysis, errors of the form factor parameters have been 
considered and inclnded in the constraint plots. The origin of \/2 comes from the 
wavefnnction of 7r° = + dd). Considering the isospin symmetry, uu = dd = 0, 

so TT^ = \/20. Using the eqnations of motion and by mnltiplying the SM hadronic 
cnrrent (H5|) by k'' and q^, see Appendix (B), we have 

Ft = (47) 

One can hnd the details of the decay rate calcnlations in Appendix (A). We hnd 
that TSM = 5.5 x 10“^^ GeV. The total decay rate of r is Tfot = 2.27 x 10“^^ GeV, 
so that BR{t~ ^ z/t- + 7 r“ + 7 r°) is 24.23% in onr calcnlations which is close to the 
experimental resnlt (25.52 ± 0.09)% [30]. Using the GVG hypothesis, it is predicted 
that BR{t- Ti-n^Ur) = (24.75 ± 0.38)% [29]. 

From the constraint t~ 7T~7r^i>r, we find that 0.07 < \Tl\ < 0.2 within the 
experimental 2a level while the theoretical nncertainty is considered within Icr. If we 
take the tensor conpling to be complex, the contour plot in Fig. [3] shows the allowed 
region of the real and imaginary components of the coupling for the measured r —>■ 
7T~7r^i>j. within the experimental 2a level and the theoretical uncertainty within lex. 
The SM expectation for the branching ratio is not allowed within the experimental 
range at the 2a level but it is allowed at higher standard deviation level. 

In the explicit leptoquark models, S'L(mLQ) = ±4TL(mLQ), one can obtain the 
constraint on Sl and Tl from r“ —)■ z/t-tt" and r“ —)■ tt^tt^z/t- at the same time. It 
is found that the limits of 0.07 < \Tl\ < 0.2 and 0.62 < \Sl\ < 1.73 are obtained 
within the experimental 2a level with the theoretical uncertainty within la. The 
allowed regions of the real and imaginary components are shown in the contour plot 
in Fig. m 

The 7171 state is produced dominantly though an intermediate vector resonance 
and is in a P wave. Therefore, the scalar terms with the couplings Sl and Sr do not 
contribute to the decay process t~ —)■ z/t-tt" as the scalar hadronic current vanishes 
because of parity. Isospin symmetry also results in Sl,r not contributing within 
decay. 

In the V ± A case, the couplings can be constrained by both r“ —)■ z/,- + 7r“ and 
r —)■ 7r“ + 7r° + z/,- decays. Gonsidering the first process the branching ratio is given 
as 


^^v±A — 5 ( 48 ) 

where the V ± A contribution is 

'■7±a = U-Vh. (49) 

From the second process, the branching ratio is given as 

+’"v'±a)^ 5 (50) 
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Figure 3: The allowed region for the real and imaginary components of the complex 
leptoquark coupling T^. The constraint on Tl is from t~ —>■ 




Re[5z,(mT)] 


Re[ri(mP] 


Figure 4: The allowed regions for the real and imaginary components of the lep¬ 
toquark running couplings SLirrir) and TL^nir) with SiirriLQ) = ±4 TiijnLQ) at 
mLQ = 1000 GeV. The constraint on Siimr) is from t~ —)■ and Ti^mr) is 

from T~ 7T~7T^Ur. 

where the V ± A contribution is 

= Vl + Vr. (51) 

If we take the couplings to be real, the contour plot in Fig. [5] shows the allowed 
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region. The allowed regions for the real and imaginary parts, if the couplings are 
taken to be complex, are shown in the contour plot in Fig. [6l 

Note, even though we consider complex couplings in the constraint equations in 
this section, we take the couplings to be real for the scattering calculations. 



Figure 5: The allowed region (green area) for the left- and right-handed couplings 
Vl and Vr- The constraints are from t~ —)• TT~Vr and t~ —)• 



Reli/J 


1.0 

0.5 

^ 0.0 
E 


- 0.5 


- 1.0 


- 2-10 1 2 
Re[F;i] 



Figure 6: The allowed regions for the left- and right-handed complex couplings Vr 
and Vr. The constraint are from t~ —and t~ —Left panel; we 
take Vr = 0 and treat Vr as a complex coupling. The allowed region is the contour 
ring. Right panel; we take = 0 and treat Vr as a complex coupling. The allowed 
region is the overlap area between the two contour rings. 
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4 Numerical analysis 

In this section the sensitivity of the neutrino cross-section scattering to the scalar and 
tensor interactions, explicit Leptoquark model, and V ± A interactions is discussed. 
We study the ratio of the total cross section, da/dxdy, and da/dt for the tan-neutrino 
to the muon-neutrino scattering. We also show the results of the total cross section, 
da/dxdy, and da/dt for the process Ur + N ^ r + X. 

4.1 Scalar and Tensor Interactions 

The ratio of the total cross section is shown in Fig.[7]while the ratio of the differential 
cross sections da/dxdy and da/dt are given in Figs. (|H1 E]). The impact of the new 
physics is clearly detectable in the ratio of the total cross section and the differential 
cross sections. The new physics effect is also observable in the total cross section, 
da/dxdy, and da/dt for the process Ur + N ^ r + X, as shown in Figs. flTUl [TTl [T^ . 

4.2 Explicit Leptoquark Model 

Here we take /ulq = 1 TeV. In Figs. flT51 HU [15]), we show the differential cross sec¬ 
tion da/dt, its ratio, and the total cross section for the particular models 5'i(mLQ) = 
±4 TL^mi^q). The impact of the new physics is clearly detectable. 

4.3 V zt A Interactions 

The ratio of the total cross section, da/dxdy, da/dt are shown in Figs. fITBl [T71 fTK]) . 
respectively. The hgures show that the effect of H ± 4 new physics is small in the 
neutrino cross section. The new physics effect is small in the total cross section, 
da/dxdy, and da/dt for the process Ur + N ^ r + X, as shown in Figs. flT^ 1201ET]) . 
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Figure 7: 5 ± T model: The ratio between the total cross section of z/,- + —)■ 
T + X to + N —>■ n + X with Scalar-Tensor couplings. The green solid 
line corresponds to the standard model predictions Sr = Sl = = 0. The 

blue dashed, red dotdashed and black dotted lines correspond to {Sr, Sl,Ti) = 
(-0.19,0.68,0.072), (1.98,0.42,-0.13), (-1.87,-1.31,0.18). 

5 Conclusion 

In this paper we discussed tests of lepton non-universal interactions through Uj. 
scattering. We adopted an effective Lagrangian description of new physics and 
considered explicit leptoquark models for our calculations. The parameters of the 
new physics were constrained by single pion and two pion r decays, T~{ki) —?• 
i^T{k 2 ) + 7r“(g) and r(p) —)■ vr“(pi) -|- vr°(p 2 ) + ^rips) , which are well measured. We 
then discussed the ratio of the total and differential cross sections for the two deep 
inelastic scattering processes + Sf ^ t + X and p^j, + N ^ + X as a probe 

of the new physics in the neutrino cross-section experiments. In the ratio of cross 
sections, the uncertainty of the parton distribution functions is expected to cancel 
out leading to precise results. In the effective Lagrangian framework we looked 
at models with scalar and tensor interactions. As an explicit realization of such 
models we considered leptoquark models where scalar and tensor couplings arise 
with relations between the couplings. Our results showed signihcant new physics 
effects, both in the total cross sections as well as in the differential distributions for 
Pr + N —)■ T + X, are allowed with the present constraints. These new physics effects 
could be observed at future proposed p^ scattering experiments. We also considered 
vector-axial vector new physics operators in our analysis. The results showed that 
the new physics effect is small in this case. 
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Figure 8: S' ± T model: The ratio between the differential cross section 
{da/dxdy) oi Vr + N —r + Xtoz/^ + iV y + X with Scalar-Tensor cou¬ 
plings. The green lines correspond to the standard model predictions Sr = 
Sr = Ti = Q. The blue, black, and red lines correspond to {Sr, Sl,Tr) = 
(-0.19,0.68,0.072), (1.98,0.42,-0.13), (-1.87,-1.31,0.18). The blue and green 
dashed lines correspond to {x,y) = ( 0 . 95,74 -|- B). The black and green dotdashed 
lines correspond to {x,y) = ( 0 . 475,(74 -|- B)/2). The red and green dotted lines 

correspond to {x, y) = { ^MiEj-m^) ^ A-B). 
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Figure 9: S' ± T model: The ratio between the differential cross section {da/dt) 
of Ur + N —>■ T + X to + N — )■ y + X in the Scalar-Tensor model. The 
green dashed, dotted and dotdashed lines correspond to the standard model pre¬ 
dictions Sr = Sr = Tr = 0 at = 30,20,10 GeV, respectively. The 
blue dashed, black dotted, and red dotdashed lines correspond to {Sr, Sr,Tr) = 
(-0.19,0.68,0.072), (1.98,0.42,-0.13), (-1.87,-1.31,0.18) at = 30,20,10 
GeV, respectively. The physical regions of the momentum transfer is taken to be 

Ql{w,ut) <Q^< Ql{w,ut). 


k' 


, '. »'> •. 

■ ■Xs \ 

••• 

\ 'Sv- 


10 20 30 40 

02 [GeV2] 


50 


17 









£^[GeV] 

Figure 10: S' ± T model: The total cross section oi Vr + N —>■ 
T + X vn the Scalar-Tensor model. The green solid line corresponds to 
the standard model predictions Sr = Sl = Tl = 0. The blue 
dashed, black dotted and red dotdashed lines correspond to {Sr,Sl,Tl) = 
(-0.19,0.68,0.072), (1.98,0.42,-0.13), (-1.87,-1.31,0.18). 



5 10 15 20 25 30 

^viGeV] 

Figure 11: S±T model: The differential cross section (da/dxdy) of Ur + N ^ r + X 

in the Scalar-Tensor model. The green lines correspond to the standard model 

predictions Sr = Sr = Tl = 0. The blue, black, and red lines correspond to 

{Sr,Sl,Tl) = (-0.19,0.68,0.072), (1.98,0.42,-0.13), (-1.87,-1.31,0.18). The 

blue and green dashed lines correspond to {x,y) = (0.95, A -f B). The black and 

green dotdashed lines correspond to {x,y) = (0.475, {A + B)/2). The red and green 

2 

dotted lines correspond to {x, y) = ( 2 M(.^-m ) ’ Si. — B). 
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Figure 12: S ±T model: The differential cross section [da/dt) of u^ + N —)■ r + X in 
the Scalar-Tensor model. The green dashed, dotted and dotdashed lines correspond 
to the standard model predictions Sr = Sl = = 0 ed = 30(left), 20(middle) 

and lO(right) GeV, respectively. The blue dashed, black dotted, and red dotdashed 
lines correspond to {Sr, Sl,Tl) = (-1.87,-1.31,0.18) at Ei, = 30(left), 20(middle) 
and 10(right) GeV, respectively. The physical regions of the momentum transfer is 
taken to be Q‘i{Wcut) < < QXiWcut)- 
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Figure 13: Leptoquark: The ratio between the differential cross section (da/dt) 
of Ur + N —)■ T + X to + N p -|- X in the Leptoquark model with 
>S'l(?^lq) = ±4 TR { mi ^ q ) at tulq = 1000 GeV. The green dashed, dotted and dot¬ 
dashed lines correspond to the standard model predictions Sr = Sl = Tl = f) aX. 
El, = 30, 20,10 GeV, respectively. The blue, black, and red lines correspond to 
{RelSiirriLQ)], Im[SL{rnLQ)], Re[TL{TnLQ)], Im[TL{mLQ)]) = (0.56, 0.60, 0.14, 0.15) 
at El, = 30, 20,10 GeV, respectively. The physical regions of the momentum transfer 
is taken to be Q'i{Wcut) < Q‘^ < Q+iWcut)- 
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Figure 14: Leptoquark: The differential cross section (da/dt) of Ur + N 
r + X in the Leptoquark model with S'i(mLQ) = ±4 TL^mi^q) at /ulq = 1000 
GeV. The green dashed, dotted and dotdashed lines correspond to the stan¬ 
dard model predictions Sr = Sl = = 0 at = 30(left), 20(middle) 

and 10(right) GeV, respectively. The blue, black, and red lines correspond to 
{Re[SLimLQ)], Im[SLirnLQ)], Re[TL{rnLQ)], Irn[TL{rnLQ)]) = (0.56, 0.60, 0.14, 0.15) 
at Ely = 30(left), 20(middle) and lO(right) GeV, respectively. The physical regions 
of the momentum transfer is taken to be Q^{Wcut) < < Q+(Wcut)- 
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Figure 15: Leptoquark: The total cross section in the Leptoquark model with 
SLimi^o) = ±4 TR^mi^Q) at /ulq = 1000 GeV. The green line corresponds to the 
standard model predictions Sr = Sr = Tr = 0. The blue line corresponds to 
{Re[SL{rnLQ)], Im[SL{rnLQ)], Re[TL{mLQ)], Im[TL{mLQ)]) = (0.56, 0.60, 0.14, 0.15). 
The physical regions of the momentum transfer is taken to be Q‘^{Wcut) < < 
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Figure 16: V±A model: The ratio between the total cross section of z/^ + iV t+X 
to + N —>■ /i + X in the V ± A model. The green line corresponds to the 
standard model predictions = Vr = 0. The red dotdashed line corresponds to 
(VL,f4) = (0.016,0.006). 
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Figure 17: V±A model: The ratio between the differential cross section (da/dxdy) 
ot Ur + N ^ T + X to i>fj_ + N ^ y + X in the V ± A model. The green lines 
correspond to the standard model predictions Vl = Vr = 0. The blue, black, and 
red lines correspond to (Vr, Vr) = (0.016,0.006). The blue and green dashed lines 
correspond to {x, y) = (0.95, A + B). The red and green dotdashed lines correspond 
to {x,y) = (0.475, {A + B)/2). The black and green dotted lines correspond (x,|/) = 
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Figure 18: V iL A model: The ratio between the differential cross section [da/dt) 
of Ur + N ^ T + X to Ufj_ + N ^ + X in the V ± A model. The green lines 
correspond to the standard model predictions 14 = Vr = 0. The blue, black, and 
red lines correspond to (14,14) = (0.016, 0.006) at = 30, 20,10 GeV, respectively. 
The physical regions of the momentum transfer is taken to be Q‘^_{Wcut) < < 



Figure 19: V E A model: The total cross section of Vr E N ^ r + X in the V E A 
model. The green line corresponds to the standard model predictions 14 = 14 = 0. 
The red dotdashed line corresponds to (14,14) = (0.016,0.006). 


22 






25 

20 

O 

o 15 

I 10 


O' 

5 10 15 20 25 30 

£v[GeV] 

Figure 20: V±A model: The differential cross section {da/dxdy) of Ur + N ^ t + X 

in the V ± A model. The green lines correspond to the standard model predictions 

Vl = Vr = 0. The blue, black, and red lines correspond to (14, Vr) = (0.016, 0.006). 

The blue and green dashed lines correspond to {x, y) = (0.95, A + B). The red and 

green dotdashed lines correspond to {x,y) = (0.475,(4 + B)/2). The black and 

2 

green dotted lines correspond to {x, y) = { 2 M{E^-m ) ’ ^ ~ B). 




Figure 21: 1/ ± 4 model: The differential cross section (da/dt) of z/,- + iV —)■ r + X 
in the V ± A model. The green lines correspond to the standard model predictions 
Vl = Vr = 0. The blue, black, and red lines correspond to (14,14;) = (0.016, 0.006) 
at Ey = 30, 20,10 GeV, respectively. The physical regions of the momentum transfer 
is taken to be Q‘i{Wcut) < Q'^ < Q+{Wcut) 
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Appendix (A) 

Here, we give details of the calculations of the process r(p) —)• i^riPs) + + 

n^{p2). In the rest frame of 7r“ and vr*^, 

P={E, P),Pi = {E,, A),P2=(^2, -Pi),P 3 = {Es, P), (52) 

k = Pi-p2= {Ei- E2, 2 Pi) ,q = Pi+P2= {Ei + E2, 0 ) , (53) 

and we define two variables, 

mu = (Pl+P2)^ = g^ 

mis = iP2+P3f- (54) 

Then, 

dr = -—^- -Xdml^dmln, (55) 

( 27 r )3 32 m 3 23 > v ; 

with 

X = - ^2 \Msm + Mrp = - ^2 + - ^2 (56) 

spin spin spin 

where Msm and Mt are given in Eqs. fl43l 1441) . and the cross terms are zero. 

By averaging the spin, we get 


XsM = iGlV^,F^{Q-^)\2{k ■p){k- Pa) - 4"(p ■ p,)], (57) 


and, 

Xt = IGGlV^^Tl ^ ^ [{k ■ qfi-p ■ ps) + 2{k ■ q){{k • p 3 ){p ■ q) 

+ {k ■ p){p3 ■ q)) - 2 q^{k ■ p){k ■ ps) + k‘^{q^{p ■ P3) - 2 {p ■ q){p3 ■ g))]- 

(58) 


All these X's can be expressed in terms of m ^2 and m 23 because 

Ei = 

E2 = 

E3 = 

E = 


ml2 

-ml 

+ m2 


2mi2 


ml2 

9 

— 

+ ml 


2mi2 


M2 . 

-m?2 

-ml 


2mi2 


m2 + m22 

1 

CO to 


2mi2 
,2 2 2 


2pi ■ p = 77123 - m2 - m3 - 2 ^ 2E3 , 


(59) 

(60) 

(61) 

(62) 

(63) 
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1^ = - h'P, 


( 64 ) 

\n\ = ^jEl-ml (65) 

with M = rur = 1.77 GeV, mi = m.,,- = 0.140 GeV, m2 = m^^o = 0.135 GeV, 

ms = m^^ = 0. 

Let’s work on the SM case first and set 


Ai = k ■p = E{Ei - E2) - 2 pi-p 

= E(^Ei — E2) + 2 E 2 E^ + m^ + m"^ — ^231 
A 2 = k-ps = Es{Ei - E2) - 2 pi ■ p 

= E^{Ei — E2) + 2E2E3 + ml + ml — ml^, 
A, = k^ = {Ei-E2y-Ap-l^ 

= {Ei-E2f-A{E!-ml), 

Ai = p-pz = EiEii - f 

= EiE^ -E‘^ + Ml 


Then, 

XsM = iGlv;i^E\Q^)[ 2 AiA 2 - 7I37I4]. 

Now, let iis integrate Xsm by within the limits 


(m23)mas. — {^2 + 



( 66 ) 

( 67 ) 

( 68 ) 

( 69 ) 

( 70 ) 

( 71 ) 


— (It'S + 




where m, — m2 — m,,, mj^ — and M — m^. One gets 


( 72 ) 


rsM — 


cose^ dg2 

96(2.)3 2 ml 



Now, we can integrate over 77142 within the limits 

Qmax ~ {A^l 2 )max = (d7 — 777i) , 

Qliin = {ml 2 )min = (^i + 7772)1 



3/2 


( 73 ) 


( 74 ) 
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Now let’s work on the tensor leptoqnark case, we set 



Bi = 

k-q = El- El, 


( 75 ) 


B 2 = 

P-P3 = EEs-'f, 


( 76 ) 


B3 = 

k-p3 = EsiEi - E 2 ) - 2 p 

■Pu 

( 77 ) 


B,= 

p - q = E{Ei + E'^), 


( 78 ) 


B,= 

k-p 3 = E{Ei - E 2 ) - 2 p- 

Pu 

( 79 ) 


B,= 

P 3 ■ q = E 3 {Ei -|- E 2 ), 


( 80 ) 


B,= 

q^ = iE, + E 2 r, 


( 81 ) 


B 8 = 

e = {E,-E 2 f-Ap-[\ 


( 82 ) 

Then, 





A't = SGlv:-,TlFUq^) r-Bj 

B 3 + 2i?i(i?3i?4 -|- BiBq) — 

2 B 1 B 1 B 3 i?8(i?7i?2 

— 2 BiBe)] 





( 83 ) 

Numerically, 

we can get Pt = 

= 3.43 X 10-12^2 QgY_ 




Appendix (B) 

In the decay process r“ —)■ z/,- + 7r“, the SM hadronic cnrrent is given in Eq. [Ml By 
mnltiplying the cnrrent by + Pu^ can End the NP scalar current given 

by 

{0\d{As - BsJ^)u\n{q)) = (84) 

If one multiplies the SM current by = p^ ~ Pu^ the scalar current will be 

(0|d(A5 - Bs'y^)u\7r{q)) = iV2f^{mu + md)Bs. (85) 


Now, by multiplying the two equations above and taking the square root, we end 
up with scalar current that is independent of the quark masses 

{0\d{As - Bs'y^)u\7r{q)) = iV2f^m^Bs. (86) 


In the process r“(p) —)■ 7i~(pi) + {p 2 ) + h'rips), the NP tensor current is given in 

Eq. 061 Here p^ = p^+Pq and P 2 = P^—Pq, where Pu and pd are the momenta of the 
up and down quarks that come from the tan decay, and (p^, —Pq) are the momenta 
of the quark-antiquark pair from the vacuum that pair up with the up and down 
quarks to form 7r° and ti~ . By multiplying the current by q^ = Pi + p 2 = Pd ~ Pu 
and using the equation of motion, in the isospin symmetry limit, one gets the form 
factor 


qi 


( 87 ) 
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If one multiplies the tensor current by —p 2 = Pd ~Pu~k ‘^Pqi the form factor 

will be given by 


Fj' —- 

{rud + mu) 


-iF 

(1-2!^) 


m —m 


rrid-mu 


( 88 ) 


Now if the tttt is dominantly coming from a vector resonance then we can expect that 
the distribution of the momenta of the quarks inside the resonance will be peaked 
around Pu = Pd- In this limit the second term in the denominator above vanishes 


as — 2^^ j = 0 . Hence,by taking the second term in the denominator small, we 
get 


Ft — — i- - -F. (89) 

[mu + md) 

Now, by multiplying the two equations above and taking the square root, the form 
factor will be independent of the quark masses 


Ft — 


-iF 

7 ?' 


(90) 
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